Field-induced 3- and 2-dimensional freezing in a quantum spin liquid. 
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Field-induced commensurate transverse magnetic ordering is observed in the Haldane-gap com- 
pound Ni(CsDi4N2)2N3(PF6) by means of neutron diffraction. Depending on the direction of applied 
field, the high-field phase is shown to be either a 3-dimensional ordered Neel state or a short-range 
ordered state with dominant 2-dimensional spin correlations. The structure of the high-field phase 
is determined, and properties of the observed quantum phase transition are discussed. 
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The one-dimensional (ID) integer-spin Heisenberg 
antifcrromagnet (AF) demonstrates unique quantum- 
mechanical behavior that is inconsistent with the conven- 
tional semi-classical model of magnetism. Due to zero- 
point quantum spin fluctuations, that destroy long-range 
order (LRO) even at T = 0, the ground state is a spin- 
singlet with short-range (exponentially decaying) spatial 
spin correlations, and is often referred to as a "quantum 
spin liquid" . The main feature of the magnetic excitation 
spectrum is the so-called Haldane energy gap pj. Quite 
remarkable is the behavior of ID and quasi-lD integer- 
spin AFs in applied magnetic fields. The effect of the 
field is to suppress zero-point fluctuations, and restore a 
gapless spectrum. The result is a quantum phase tran- 
sition at a certain critical field H c , to a Neel- like state 
with long-range order that may be characterized as a 
"spin solid". Thus, in an unusual twist, a uniform field 
induces staggered magnetization in quantum-disordered 
spin chains. 

The magnetization process of a S = 1 ID AF is 
now rather well understood theoretically pj-f/j. For the 
isotropic case the problem was shown to be equivalent 
to Bose condensation in one dimension j|]g]. Many the- 
oretical results were confirmed in experimental studies 
of real quasi-lD compounds. For a long time however, 
the actual phase transition remained inaccessible for ex- 
perimental investigation. In some materials (such as 
Y2BaNiC>5, for example ||), the value of H c is pro- 
hibitively high. In other compounds {e.g., NENP pPf), 
the transition does not occur due to certain structural 
features, and is instead replaced by a broad cross-over 
phenomenon ]ll|-|l3[|. A true phase transition at H c has 
been experimentally observed only recently, in the quasi- 
lD S = 1 AF materials NDMAP @-[l|] and NDMAZ 
p7[ . Specific heat, magnetization and ESR studies have 
provided a comprehensive picture of the H — T phase dia- 
gram IliMlq] , a refined version of which is shown in Fig. [jL 
Nevertheless, to date, the nature of the high-field phase 
remained undetermined, and no direct evidence of stag- 
gered LRO has been obtained. In particular, preliminary 



neutron scattering experiments of Rcf. JL8| failed to de- 
tect any magnetic Bragg peaks above H c . In the present 
work we report a high-field magnetic neutron diffraction 
study of NDMAP. We find that the high-field phase is 
characterized by a commensurate ordering of spin com- 
ponents perpendicular to the field direction. Surprisingly, 
the nature of the high-field phase is extremely sensitive 
to direction of applied field: depending on the experi- 
mental geometry, for H > H c the system either develops 
true 3D LRO, or static quasi- 2-dimensional short range 
order. 
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FIG. 1. Field-temperature phase diagram of NDMAP. 
A quantum-disordered (QD) phase with a gap is seen in 
low fields. The high-field phase is characterized by static 
long-range order (3D-LRO) or quasi-2D short-range order 
(2D-SRO). Open symbols from specific heat measurements, 
as in Ref. |I5| . Solid symbols from neutron diffraction. 

Magnetic properties of NDMAP are due to S = 1 Ni 2+ 
spins that are assembled into AF chains running along 
the c axis of the orthorhombic Pnmn crystal structure 
(a = 18.05 A, b = 8.71 A, c = 6.14 A). Haldane gap exci- 
tations and magnetic interactions in this system were pre- 
viously investigated by means of zero-field inelastic neu- 
tron scattering |ljj. In-chain AF interactions are domi- 
nant and correspond to a Heisenberg exchange constant 



J = 2.8 meV. The triplet of Haldane excitations is split 
by easy- (a, b) plane anisotropy, and the gap energies are 
A z = 1.9 meV and A_l « 0.47 meV, for excitations po- 
larized along and perpendicular to the anisotropy axis, 
respectively. The coupling between chains along the b 
axis is rather weak, J y = 2 • 1CP 3 meV, and barely de- 
tectable along a, \J X \ < 5 • 10~ 4 meV. 

The present neutron diffraction experiments were per- 
formed at the SPINS cold-neutron 3-axis spectrometer at 
the National Institute of Standards and Technology Cen- 
ter for Neutron Research. A 150 mg sample was mounted 
with either the (0, k,l), (h,h,l) or (h,k,k) reciprocal- 
space planes in the scattering plane of the spectrometer. 
Most measurements were performed using ( 58 Niguide) — 
80' - 80' - 240' or ( 58 Niguide) - 40' - 40' - 240' colli- 
mations, with a Be-filter positioned in front of the sam- 
ple to eliminate higher-order beam contamination, and 
Ef = 5 meV or Ef = 3 meV fixed final energy neutrons. 
Sample environment in all cases was a pumped He-3 cryo- 
genic insert and a 9 T or 7 T superconducting magnet, 
with the field applied vertically, i.e., perpendicular to the 
scattering plane. 




HI [0,1,-1] 


00 


q=(<U1 




SP$rO 







h (r.l.u.) 



0.50 

/(r.l.u.) 



H|[l, 


,0] (c) 


q-(A,-0 .5,1.5) 


5.J.S, 


u 


/ 


?4 


/ 




/ 


\ 


J 


\ 





H II [1.0.0] 


(d) 


q-(0,*,l-5) 




°A 








a^y 


J?o 








H || [1,0.0] 
. q=(0,-0.5,f) 



(e) 




-0.2 0.0 0.2 -0.55 -0.50 -0.45 1.40 1.45 1.50 1.55 1.60 

h (r.l.u.) A (r.l.u.) /(r.l.u.) 



FIG. 2. Typical elastic scans (open symbols), showing 
magnetic Bragg scattering in NDMAP at T = 0.25 K, 
H = 9 T||[0, -1,1] (a,b) and T = 0.375 K, H = 7 T]|[l, 0, 0] 
(c-e). The solid lines and symbols represent experimental 
wave vector resolution. Dashed lines as described in the text. 

True 3D commensurate AF LRO was observed in a 
crystal mounted in the (h, k, k) zone, with the magnetic 
field applied along the [0,-1, 1] direction. At T = 0.25 K 
the transition occurs at H c = 4.8 T, and manifests 
itself in the appearance of new Bragg reflections, at 
(h, 2 i 2 ) (h-even integer; n- integer) reciprocal- 



space positions. No additional peaks were found at 
odd-h reciprocal-space points. Figures ||(a) and ||(b) 
show typical scans through the (0,0.5,0.5) peak mea- 
sured at H — 9 T applied field (symbols). The peaks are 
resolution-limited along both the (1, 0, 0) and (0, 1, 1) di- 
rections. Wave vector resolution (Fig. ||, solid lines) was 
either calculated or directly measured in our experiment. 
The multiple peaks seen in the rocking scan in Fig. 0(a) 
result from sample mosaic. Geometrical limitations im- 
posed by the use of the bulky 9 T magnet prevented 
us from measuring the width of the peak in the vertical 
[0, 1,-1] direction. However, the substantial values of the 
b and c axis coupling constants, as well as the resolution- 
limited peak widths along a, where magnetic interactions 
are weakest, suggest that the observed reflections indeed 
have zero intrinsic width in all three directions, and rep- 
resent true 3D long-range order. This is also consistent 
with the measured absolute values of magnetic intensi- 
ties, as discussed below. 

To determine the spin arrangement in the high- 
field ordered phase, intensities of (0,0.5,0.5) (2,0.5,0.5) 
(4,0.5,0.5) (2,1.5,1.5) and (4,1.5,1.5) magnetic Bragg 
peaks in the (/i, k, k) plane were measured in rocking 
scans at T — 0.25 and H = 9 T. Magnetic intensities 
were brought to the absolute scale by comparing them to 
nuclear intensities measured in the same configuration at 
room temperature (the exact low-temperature structure 
of NDMAP has not been determined to date). The data 
were analyzed using a simple collinear model. As the en- 
ergy scale of in-chain AF interactions is much larger than 
that defined by the applied magnetic field ( J> g/iB-ff), 
the ordered moment is expected to be perpendicular to 
the field direction. The (a, 6)-easy plane anisotropy en- 
sures that it should also be perpendicular to the chain- 
axis. For a field applied along [0,-1, 1] the spins are thus 
along (1,0,0), as confirmed by recent ESR experiments 
||15i . This allows only one magnetic structure that would 
be consistent with the observed selection rule for mag- 
netic Bragg peaks. Varying only the sublattice magneti- 
zation m, we obtained a good fit to the data (x 2 = 1.6) 
with to = 1.13(5) /ib- This value is to be compared to 
to = 2 /!b for a fully saturated S = 1 antiferromagnet. 

The measured field dependence of the (0, 0.5, 0.5) mag- 
netic Bragg peak intensity at T = 0.25 K is shown in 
Fig. [|(a) (solid circles), and is consistent with our ex- 
pectations for an easy-plane Haldane system at T — > 0. 
The data were analyzed assuming a power-law field de- 
pendencies of the Bragg intensity and ordered moment: 
1(H) ex (H - H c ) 213 . A fit of this equation to the 
data (Fig. |(a), solid line) yields H c = 4.81(1) T and 
j3 = 0.207(5). Since the transition at H c occurs through 
a softening of the lowest-energy Haldane gap excitation 
B, H c is determined by the Haldane gap energies. It 
is straightforward to show that for a magnetic field ap- 
plied at angle a to the anisotropy axis c is given by: 



H B H C = A z Aj_/Jg 2 A 2 z cos 2 a + g 2 ± A z A ± sin 2 a. In this 
formula g z = 2.1 and g±_ = 2.17 are components of the 
Ni 2+ gyromagnetic tensor along and perpendicular to 
the anisotropy axis, respectively [Q. For H ||[0, — 1, 1], 
a = 54.8° and this equation gives H c — 5.4 T, in rea- 
sonable agreement with our experimental result. In the 
case of broken axial symmetry, for a single chain in a 
magnetic field at T = 0, the transition is expected to 
fall in the 2D Ising universality class, with magnetic field 
taking the role of effective temperature B. The order 
parameter critical exponent for this model is (3 — 0.125. 
For NDMAP however, inter-chain interactions along the 
b axis can be considered substantial, and 3D Ising be- 
havior, with /3 « 0.31, may be expected. The measured 
critical exponent falls in between these two values and is 
characteristic of a dimensional crossover regime. 

The temperature dependence of the (0, 0.5, 0.5) mag- 
netic Bragg intensity measured at H = 9 T is shown in 
Fig. 0(b) (solid circles) . An analysis of the data assum- 
ing a power law (solid line) yields T c = 1.7(1) K. The 
critical exponent was found to be indistinguishable from 
the mean field value (3 = 0.5, indicating that at T w 1.5 K 
the true critical region is too narrow to be investigated 
in the present experiment. 

A striking result of this work is that for a mag- 
netic field applied perpendicular to the chain-axis, the 
high-field phase is no longer a 3D-ordered state, but 
has predominantly 2-dimensional spin correlations. For 
the crystal mounted in the (0, fc, I) zone, at T = 
0.38 K and JET||[1,0,0], magnetic scattering was de- 
tected at (0, 2 " 2 +1 , 2 "i +1 ) (n, m-integer) reciprocal-space 
points above H c = 6 T. k- and Z-scans through the 
(0,-0.5,1.5) reflections contain well-defined peaks, as 
shown in Figs. g(d) and g(e) (symbols). Scans perpendic- 
ular to the (&, c) plane however, reveal that the scattering 
is concentrated in Bragg rods along the (1, 0, 0) direction, 
rather than Bragg peaks, as in the -ff||[0, —1, 1] case. Fig- 
ure. 0(c) shows an /i-scan through the (0, —0.5, 1.5) posi- 
tion, that could only be performed in a limited range 
due to geometrical constraints. The measured inten- 
sity is independent of h. Note that the measured ver- 
tical resolution (Fig. u(c), solid line), is sufficient to 
observe a well-defined Bragg peak. The rod nature of 
magnetic scattering implies that static spin correlations 
along the a axis are absent in the system, and that mag- 
netic ordering occurs only within individual (6, c) planes. 
This is consistent with the a-axis being the direction of 
weakest magnetic interactions. Similar behavior was ob- 
served for a magnetic field applied along the [1,-1,0] 
direction, where magnetic scattering was observed above 
H c = 5.8(2) T, at T = 0.25 K, and also found to be 
concentrated in Bragg rods parallel to the a axis. In 
this geometry, the rods cross the (h, h, I) scattering plane 
±,2itt!,^ti) ^ m _ integer ) positions. The 



netic diffraction intensity at q = (0.5, 0.5, 0.5) are shown 
in Fig. (open symbols). A well-defined transition is 
observed in agreement with specific heat measurements 
(Fig. [l], diamonds). 
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FIG. 3. Measured field (a) and temperature (b) depen- 
dence of the ordered moment in NDMAP for a magnetic field 
applied in the (b, c) plane. Solid lines are power-law fits to 
the data, as described in the text. Inset: A schematic view of 
the spin arrangement in the high-field phase. 

The apparent Bragg rod intensity for H\\ [1,0,0] and 
H 1 1 [1,-1,0] is substantially smaller than the intensity 
of magnetic Bragg peaks seen for iJ||[0, 1, —1] (Fig. ||). 
This is due to that fact that in the 2D-ordered phase 
magnetic intensity is spread out through an entire Bril- 
louin zone along the a axis, rather than being concen- 
trated at h = 0. However, in the 2D phase, within each 
(6, c) plane, the actual ordered moment is similar to that 
in the 3D ordered state. The 2D ordered moment for 
i?||[l,0,0] and ff||[l,— 1,0] was deduced from an anal- 
ysis of measured Bragg rod intensities, paying close at- 
tention to resolution effects associated with vertical an- 
gular acceptance of the spectrometer. A series of rocking 
scans across the rods were collected for H = 7 T| | [1, 0, 0] , 
T = 0.38 K and H = 9 T||[l,-1,0], T = 0.25 K, re- 
spectively. The data were normalized by the measured 
nuclear Bragg peak intensities. The 2-dimensional spin 
structure was assumed to be collinear, with spins per- 
pendicular to both the chain axis and H in each case. 
Nearest-neighbor spins were assumed to be aligned an- 
tiparallel along both the b and c axes. This simple model 
was found to agree very well with the available data. The 
ordered moment was determined to be m = 0.66(4) \xb 
and m = 1.2(1) fi B for H = 7 T||[l, 0,0], T = 0.38 K 
and H = 9 T||[l,-1,0], T = 0.25 K, respectively. The 
latter value agrees very well with the 3D ordered moment 
measured at H ||[0, — 1, 1] at the same value of field and 
temperature. 

A careful analysis of the scans across the Bragg rods 
for W| | [1,0,0] reveals a finite correlation length even 



within the quasi-2D ordered (b, c) planes. Along the di- 
rection of strongest coupling, i.e., along the chains, the 
magnetic peaks are resolution-limited. This can be seen 
from Fig. @(e), where the solid line represents experi- 
mental /-resolution. However, along the b axis, where 
magnetic interactions are weaker, the scattering rods are 
visibly broader than experimental resolution (Fig. @(d), 
solid line). Fitting the fc-scan to a Lorenzian profile con- 
voluted with the resolution function (Fig. ^|(d), dashed 
line) allows us to extract the intrinsic width of the 
rod, that corresponds to a real-space correlation length 
£b = 100(20) ha 106. The high-field phase in this case 
is thus characterized as a short-range ordered state with 
almost-perfect correlations along the chains and short- 
range correlations in the transverse direction. 

This type of "spin freezing" in NDMAP is very simi- 
lar to that recently found in the S — 1/2 quasi-lD AF 
SrCuC>2 |Q. For a 2-dimensional Hcisenberg system, 
the ground state is disordered and gapless for half-integer 
and odd-integer spin due to "instanton" topological de- 
fects [gl[ . If magnetic coupling along the third direction 
is sufficiently weak, as it is in NDMAP, it will not be able 
to restore LRO in the system due to pinning of such de- 
fects within each plane. The similar behavior of two very 
different systems, NDMAP and SrCuC>2, suggests that 
anisotropic spin freezing may be a rather general feature 
of quasi- low dimensional antiferromagnets. 

What makes the case of NDMAP special is that both 
true 3D LRO and spin freezing can occur in the same 
sample, depending on the direction of applied field. This 
may result from the anisotropic nature of magnetic in- 
teractions along the a axis, that, due to the very long 
lattice spacing, are expected to be largely dipolar. The 
strength and even sign of such coupling depends on the 
orientation of the ordered moment, which, in turn, is de- 
fined by the direction of applied field. Changing the field 
orientation for NDMAP is thus a way to tune those mag- 
netic interactions ultimately responsible for 3D LRO or 
spin-freezing behavior. 

For the basic physics of Haldane spin chains, a very 
important experimental result is the observation of a 
commensurate ordering of transverse (relative to H ) spin 
components above H c . Such staggered LRO in a quasi- 
lD material is a direct consequence of a divergence in 
the transverse spin correlation function for an isolated 
chain in an external field at q = n. Theory predicts 
that above H c a divergence also exists in the longitudi- 
nal correlator, but at a field-dependent incommensurate 
wave vector P,H). No incommensurate magnetic peaks 
were found in our experiments on NDMAP. Of course, a 
divergent susceptibility in a single chain does not guaran- 
tee the appearance of LRO at the same wave vector in a 
3D material. It is also conceivable that incommensurate 
elastic scattering in NDMAP appears at a different wave 
vector transfer perpendicular to the chain axis, due to 
anisotropy of inter-chain interactions, and thus escapes 



detection in the present study. Further work will be re- 
quired to fully resolve this issue. 

In summary, our results provide direct experimental 
evidence of field-induced commensurate AF order in a 
Haldane-gap system, and are in good agreement with 
theoretical expectations. Specifics of inter-chain inter- 
actions in NDMAP result in a rich phase diagram with 
the high-field phase being either a true 3D ordered state 
or a highly anisotropic short-range ordered state. 
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